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The High Energy Stereoscopic System (H.E.S.S.) is an array of five Imaging Atmospheric 
Cherenkov Telescopes (lACTs) designed to detect cosmogenic gamma-rays with very high ener¬ 
gies. Originally consisting of just four identical lACTs (CTl-4) with an effective mirror diameter 
of 12 m each, it was expanded with a fifth lACT (CT5) with a mirror diameter of 28 m in 2012. 
Being the largest lACT worldwide, CT5 allows to lower the energy threshold of H.E.S.S., mak¬ 
ing the array sensitive at energies where space-based detectors run out of statistics. Events can be 
analysed either monoscopically (i.e. using only information of CT5) or stereoscopically (requir¬ 
ing at least two triggered telescopes per event). To achieve a good performance, a sophisticated 
event reconstruction and analysis framework is indispensable. This is particularly important for 
H.E.S.S. since it is now the first lACT array that consists of different telescope types. An advanced 
reconstruction method is based on a semi-analytical model of electromagnetic particle showers 
in the atmosphere (“model analysis”). The properties of the primary particle are reconstructed 
by comparing the image recorded by each triggered telescope with the Cherenkov emission from 
the shower model using a log-likelihood maximisation. Due to its performance, this method has 
become one of the standard analysis techniques applied to CTl-4 data. Now it has been modified 
for use with the five-telescope array. We present the adapted model analysis and its performance 
in both monoscopic and stereoscopic analysis mode. 
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1. Introduction 

Within the last decade, the field of ground-based gamma-ray astronomy has evolved from an 
experimental niche with only few detected sources to a mature branch of astrophysics. The High 
Energy Stereoscopic System (H.E.S.S.) has been playing a key role in this transition process. It 
was inaugurated in 2004 with four 12m-size Imaging Atmospheric Cherenkov Telescopes (lACTs) 
located in Namibia, marking the beginning of H.E.S.S. I. The greater part of the currently known 
very high energy (VHE, E > lOOGeV) gamma-ray sources were detected with these four lACTs 
(named CTl-4). In 2012, a fifth telescope (CT5) was added in the centre of the array, initiating 
H.E.S.S. II. With its large mirror with a size of 24m x 32m, it is not only able to detect elec¬ 
tromagnetic showers which are too faint to be seen by CTl-4, but also makes H.E.S.S. the only 
lACT array that consists of different telescope types and sizes. A capable event reconstruction and 
analysis framework is hence necessary to live up to the success of H.E.S.S. I. 

In this contribution, we present the successful adaptation of the model analysis [1] to the 
full five-telescope H.E.S.S. array. The focus is set on the illustration and comparison of different 
reconstruction and analysis modes. After describing the revised model analysis in Section 2, three 
different analysis modes of H.E.S.S. II data are tested on the Crab Nebula in Section 3. After that, 
the performance is evaluated in Section 4. 


2. The Model Analysis 


In 2009, [1] introduced an advanced event reconstruction method for lACTs^ Its founda¬ 
tion is a semi-analytical model that describes the average development of electromagnetic particle 
showers in the Earth’s atmosphere as a function of the primary energy of a gamma-ray {E) and its 
depth of first interaction (T). Building up on this model, the amount of Cherenkov emission of a 
shower in the camera of an lACT is analytically calculated. Templates are generated and stored for 
different zenith angles, impact distances^ R, E, and T. Eor the event reconstruction, the intensities 
as measured by the camera of an lACT are compared to the shower templates by calculating a 
likelihood for each pixel: 
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( 2 . 1 ) 


where s denotes the measured signal intensity and /i the expectation value from the model tem¬ 
plates, the latter being interpolated between the grid points of the templates. Gp corresponds to 
the pedestal width, Gy to the width of the single photoelectron peak, and Gc accounts for possible 
calibration uncertainties. To compare measured and modelled intensities for a set of pixels (e.g. all 
pixels of a telescope or the whole array), a log-likelihood function is defined as 


luTset — ^ InT/— ^ 2 X lnP(^/|/i/, <7p, <7y, Cc) . (2.2) 

pixel i pixel i 

^This article is the baseline reference of this section. 

^The impact distance corresponds to the distance of the projected shower axis to the telescope. 
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Figure 1: Sketch that illustrates the different reconstruction and analysis modes (see text for further expla¬ 
nation). 


To reconstruct the parameters of a primary particle (assuming that it is a gamma-ray), the log- 
likelihood is optimised using a Levenberg-Marquardt algorithm (see [2] and [3]). The result of 
the optimisation is an estimate of the arrival direction Dir, E, T, and R. Additionally the uncer¬ 
tainty of each parameter is assessed by calculating the second derivatives of InL. As an example, 
AD/r = 1 / \/d^\nL/dDir^ provides an event-wise estimate on how well the direction was recon¬ 
structed. 

During H.E.S.S. phase I, at least two telescopes of the array were required to be triggered by 
a shower for the readout of an event, and thus no single-telescope reconstruction was performed. 
Since the inauguration of CT5, H.E.S.S. reads out both stereoscopic events from CTl-5 and mono- 
scopic events from CT5. The reconstruction of monoscopic events from such a large telescope 
allows to significantly lower the energy threshold. Eig. 1 shows the different reconstruction and 
analysis modes for the H.E.S.S. II model analysis. The stereoscopic events where CT5 triggered 
are reconstructed both monoscopically and stereoscopically. In addition to the Mono and Stereo 
analysis modes, we introduce a Combined analysis mode which makes use of both reconstruction 
approaches. In case an event was fitted both monoscopically and stereoscopically, the respective 
AD/r values are used to decide which result to use for the Combined mode. 

To separate cosmic gamma-rays from the excessive background which is mainly caused by par¬ 
ticle showers of charged cosmic rays, the Shower Goodness compares the optimised log-likelihood 
of the shower region with its mean expectation value (see [1] for the calculation of the latter): 


Gsg = 


L- [lnL(.,|M,) - (InL) |^,] 
\/2 ■ NdF 


(2.3) 


where NdE denotes the number of degrees of freedom and InL is calculated as defined in Eq. 2.2. 
Gsg is designed to behave like a normal variable. Eor stereoscopic reconstruction, the Shower 
Goodness values of the individual telescopes are averaged to obtain one overall Shower Goodness. 

To reject events that are compatible with night-sky background (NSB) noise, the MSB Good¬ 
ness is defined as 


Gnsb = 


L-[lnL(.,-|0)-(lnL) |o] 


(2.4) 


\/2l^ 

The main purpose of this variable is to adjust the energy threshold of the analysis, which is needed 
to limit systematic uncertainties. 
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Mono Stereo Combined 



[-4,0.6] [-4,0.9] [-4,0.9] 

>32 >28 >32 


AD/r 


T [-1.1,1.3] [-1.1,3.4] [-1.1,3.4] 

^Dir < 0.3° < 0.2° < 0.3° 

t?2(deg2) <0.015 <0.006 <0.015 


Table 1: Acceptance ranges of the respective Standard cut configurations for the three different analysis 
modes. The respective cuts were optimised for point-like sources. 

For each analysis mode (Mono, Stereo, and Combined) a Standard cut configuration was 
defined (see Table 1). The event selection cuts, except for the cut on Gnsb, were optimised for 
maximum significance on a weak point source with a power-law spectrum observed at a zenith 
angle of 18°. The photon index was set to 3.0 for optimising the Mono Standard cuts, and 2.5 for 
Stereo Standard. Additional cut configurations were defined to maxinfise the discovery potential 
for sources with other spectral indices and/or cutoffs. Systematic uncertainties related to the event 
selection cuts were included in consideration, with a minor effect on the final choice of the analysis 
cuts. The cuts for the Combined mode have not been optinfised yet but correspond to a compromise 
between the ones of the Mono and Stereo modes. 

3. Application to the Crab Nebula 

To test the three H.E.S.S. II model analysis modes, 7.47 h of good-quality observations of the 
Crab Nebula are used. After dead-time correction, this corresponds to 7.23 h {Mono and Combined 
mode) and 6.79h {Stereo mode), respectively. The data were taken in October and November 2014 
at zenith angles between 45° and 55°. The source was observed in wobble mode with offsets of 
zb0.5° in both right ascension and declination. 

The significance sky map as obtained with the ring background method [4] in Mono mode 
is shown in the left panel of Fig. 2. It shows a radially symmetric point-like source (detailed 
comparisons of the simulated and measured angular resolution still have to be performed though) 
that is centered on the position of the Crab Nebula. The map is very well normalised (see right 
panel of Fig. 2), indicating that background subtraction errors are well under control. 

Fig. 3 shows the angular distributions and event statistics of gamma-like events for all three 
analysis modes. For each mode, the source and background distributions are in good agreement at 
large values. The Combined mode yields the highest excess rate and significance whereas the 
Stereo mode performs best in terms of S/B ratio. 

4. Performance 

In the following, the performance of the three modes for low-zenith observations is evaluated 
using Monte-Carlo (MC) simulations. The simulations correspond to a point source at wobble 
offset 0.5°, zenith angle 18°, and azimuth angle 180° (i.e. southward pointing). The results of this 
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Figure 2: Left: significance map of the region around the Crab Nebula, obtained using the monoscopic 
analysis mode. Right: distribution of significances for the whole map (dark blue) and the part outside the 
target exclusion region (red histogram). The latter was fitted with a Gaussian (solid red line; see fit results 
on the top right). Only significances up to 40cr are shown for visibility reasons. 

study should not be directly compared to the ones of the previous section since the Crab Nebula 
can just be observed at zenith angles above 45° with H.E.S.S.. 

The post-cut effective area of the H.E.S.S. II model analysis for the respective Standard cuts 
of the three analysis modes is shown in Eig. 4. The effective area of CTl-4 Stereo (H.E.S.S. I) is 
drawn for comparison. Whilst being almost identical at E > 300GeV, the curves of CTl-4 Stereo, 
CTl-5 Stereo, and CTl-5 Combined diverge below. By construction, the Combined mode joins 
the low energy threshold of the Mono analysis with the high detection area of the Stereo mode at 
higher energies. It thus provides the best energy coverage. 

To assess the overall performance of the revised model analysis, the differential sensitivity was 
calculated. Here it is defined as the minimum source strength to obtain a 5(7 detection after 50 h 
of observation time at a given energy range. Eor the calculation of the significance, the simplified 
method Ny/ y^A^bkg used. Background events were taken from actual observations with zenith 

angles between 12° and 22° to closely match the simulations. The differential sensitivities for the 
different modes are shown in Eig. 5. Each of the curves was calculated using five bins per decade. 
A minimum of 10 excess counts as well as a minimum S/B ratio of 0.05 was required for each bin. 
As already expected from the effective area curves, the Combined mode always performs at least 
as good as the Mono one. At medium and high energies it is however currently not as sensitive as 
the Stereo mode, which can partly be attributed to the fact that its cuts are not yet optimised. The 
larger cut for instance directly translates into an increased background level. 

5. Conclusions 

We presented the successful adaptation of the advanced model analysis [1] to H.E.S.S. phase 
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Figure 3: Angular distribution of 7 -like events from the Crab Nebula for the monoscopic {top), stereoscopic 
{middle), and combined {bottom) analysis mode, respectively. In each panel, the distribution around the 
source is shown in dark blue, and the one from the background regions (scaled with a) in red. The cut is 
illustrated with a dashed green line, and the statistics are given on the top right. 


II. In addition to the Mono and Stereo analysis modes, we introduced the first combined analysis 
of both monoscopic and stereoscopic events for an lACT instrument. This method provides the 
best energy coverage and makes use of all the information recorded by the five-telescope H.E.S.S. 
array. 

The analysis modes were successfully tested on the Crab Nebula, and all provide trustworthy 
results on this standard source. For complementation, the effective area and differential sensitivity 
were evaluated for sources with low zenith angles using MC simulations and background events 
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Figure 4: Effective area (after selection cuts) of the model analysis for the different analysis modes presented 
here. The curves were computed using simulations with a zenith angle of 18° and an azimuth angle of 180°. 



Figure 5: Differential sensitivity of the model analysis for the three analysis modes of H.E.S.S. II data, as 
well as for H.E.S.S. I for comparison. All curves correspond to point-like sources. The reference spectrum 
of the Crab Nebula used to draw the dashed black lines was taken from [5]. 
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from actual data. Though some optimisation still has to be carried out for the Combined mode, 
all results shown here prove that the updated model analysis for H.E.S.S. II continues to be a very 
sensitive and advanced method in the field of ground-based gamma-ray astronomy. 
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